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ABSTRACT 

We study the outflow activity, photometric variability and morphology of 
three very young stellar objects in the Serpens NW star forming region: 00 Ser- 
pentis, EC 37 (V370 Ser) and EC 53 (V371 Ser). High spatial resolution 
Keck/NIRC2 laser guide star adaptive optics images obtained in 2007 and 2009 
in broad-band K and in a narrow-band fllter centered on the 1-0 S(l) emission 
line of H2 allow us to identify the outflows from all three objects. We also present 
new, seeing-limited data on the photometric evolution of the 00 Ser reflection 
nebula and re-analyze previously published data. We flnd that 00 Ser declined 
in brightness from its outburst peak in 1995 to about 2003, but that this decline 
has recently stopped and actually reversed itself in some areas of the reflection 
nebula. The morphology and proper motions of the shock fronts MHO 2218 near 
EC 37 suggest that they all originate in EC 37 and that this is an outflow seen 
nearly along its axis. We identify an H2 jet emerging from the cometary nebula 
EC 53. The star illuminating EC 53 is periodically variable with a period of 543 
days, and has a close-by, non-variable companion at a projected distance of 92 
AU. We argue that the periodic variability is the result of accretion instabilities 
triggered by another very close, not directly observable, binary companion and 
that EC 53 can be understood in the model of a multiple system developing into 
a hierarchical conflguration. 

Subject headings: binaries: close — ISM: jets and outflows — stars: formation 
— stars: low-mass — stars: variables: T Tauri, Herbig Ae/Be 



1. INTRODUCTION 

The Serpens star forming region is a compact region of low-rn ass star formation. At 



a distance of 311 pc (Ide Lara. Chavarria-K. fc Lopez-Molina I Il99ll ) it is still close enough 
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for detailed, spatially resolved studies of its individual young stellar objects (YSOs) and has 
therefore attracted considerable attention. The Serpens NW region is part of the Serpens 
molecular cloud and is the site of much of the most recent star forn iation activity in this 
cloud, as indicated by multiple sour ces of CO outflow and H2 jets (IHodapp et al.l Il996al : 
Herbst. Beckwith. fc Robbertdll997l ) in an area of less than 0.2 pc in projected diameter. 
After two decades of monitoring this region in the K band, it is clear that two objects 
in Serpens NW are particularly noteworthy in the context of establishing the relationship 
between outflows, accretion instabilities, and multiplicity: Serpens DEOS = 00 Ser, and 
EC 53 = V371 Ser. 

The S erpens NW molecular core contains the unique eruptive variable 00 Ser dis- 
covered by Hodapp et al. (jl996a ) and originally named the "Deeply Embedded Outburst 
Star (DEOS)". This object has an outburst amplitude similar to the larger and more lumi- 
nous FU Orionis (FUorl-c lass of o b jects and an outburst duration longer than that of EX 
Lupi (EXor) objects. (See iHerbig I ( 1l977l ) for the definition of these object classes.) It is so 
deeply embedded as to be inaccessib l e at optical wavelen gths and shows an infrared spectrum 
( jHodapp et al.lll996at lHodapplll999l : iKospal et al.ll2007l ) much more dominated by dust than 
either of those established classes of young eruptive variables. Based on its as sociation with 



a bipol ar nebula, and its detection at sub-mm and millimet er wavelengths. iHodapp et al. 



( Il996a[ ) have classified this object as Class I in th e schem e of lAdams. Lada. fc Shu I ( 119871 ). 
This conclusion was corroborated by iKospal et al.l ( 120071 ) who found, on the basis of archival 
Spitzer data, that 00 Ser is the dominant source in Se r pens N W at mid-infrared wavelengths. 
This point was further strengthened by lEvans et al.l (120091 ) who observed the Serpens star 
forming region with Spitzer as part of the "Clouds to Disks" (C2D) project. Among the 
sources in Serpens NW classified in their survey, 00 Ser has the lowest bolometric temper- 
ature (130 K) and highest luminosity (10 Lq). 00 Ser consists of an extended reflection 
nebula and a more compact central flux maximum. The latter appears slightly extended on 
our adaptive optics images with a cometary morphology opening to the west and therefore 
represents the most intense inner parts of the reflection nebula. At near-infrared wave- 
lengths, we probably do not see any direct, unscattered light from the embedded star. We 
therefore use the variable star designation 00 Ser to describe the variable reflection nebula 
in its entirety. 

The variable object V3 71 Ser, originally name d EC 53 in the IR imaging survey of the 
Serpens molecular cloud by lEiroa fc Casalil (119921 ) . and the irregularly variable star V370 
Ser (EC 37) have also attracted attention. Both these sources are probably among the 
more evolved YSOs in the Serpens NW region, since they are quite bright at near-infrared 
wavelengths and do not dominate at mid-IR wavelengths. V371 Ser (EC 53) is a "cometary" 
nebula, i.e., one visible lobe of a bipolar structure. Both objects are associated with molecular 
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hydrogen shock fronts indicating ongoing outflow and therefore accretion activity. 

One goal of this st udy was to deterni ine the source of the faint H2 v=l-0 S(l) emission 



labeled MHO 2218 by iDavis et all f l2010h in the area of EC 37 and 00 Ser, to determine 
which of these two objects is the driving source of this shock excited emission. In a complex, 
dense star-forming region such a Serpens NW, it is not trivial to properly associate H2 shock 
fronts with their driving sources, yet such proper identification is essential for determining the 
evolutionary state of the YSOs. Another purpose was to study the bipolar nebula associated 
with the eruptive variable object 00 Ser and the changes it underwent since its last observed 
eruption in 1995. The cometary nebula EC 53 was fortuitously included in the field of view 
and the study of this object returned very interesting results. 

We have observed the Serpens NW region with the Keck II laser adaptive optics system 
both in broad band Kg and in H2 line emission on two occasions separated by about two 
years. The observations are listed in Table 1, described in detail in section 2.1, and the 
results are discussed in section 3.1. To further investigate the temporal evolution of the light 
distribution in the 00 Ser bipolar nebula, older, seeing-limited data, primarily from the UH 
2.2 m telescope, but also recently obtained data from the Infrared Imaging System (IRIS), 
UKIRT, and the CFHT were analyzed in addition to the Keck data. These observations are 
described in Section 2.2, are listed in Table 2, and discussed in Sections 3.2 and 3.3. 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. Keck Adaptive Optics Imaging 

The adaptive optics data reported here were ob tained at the Keck II t elescope, using 



the laser guide star adaptive optics (LGSAO) system ( IWizinowich et al.ll2006l ) in conjunction 
with the Near-Infrared Camera 2 (NIRC2) camera. The dates, details of the observations, 
and the achieved FWHM image quality of the final images are summarized in Table 1. In 
addition to the object images, we took separate sky frames at positions that appeared free 
of bright stars and extended emission on deep, seeing-limited images taken earlier. The 
dithered sky frames were also adaptive optics corrected to allow the optimal removal of any 
faint stars by computing the median of a set of such sky frames. Flatfield correction was done 
using differential (lights-on minus lights-off) dome fiats. The Serpens NW region is a dense 
nearby molecular cloud and consequently suffers from a lack of available optically visible tip- 
tilt guide stars. Within the constraints of the Keck LGSAO system, the only available guide 
star is EC 92 (2MASSJ18295360+0117017) with an R magnitude of 16.6. Even with this 
tip-tilt guide star, 00 Ser can only be reached when using the wide field (40" FOV) mode 
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of NIRC2, and only with a rotator position of Ri279° to place the tip-tilt star in the extreme 
corner of the guider field. These circumstances explain the otherwise strange choice of field 
orientation and object placement in our images. It is obvious that the object placement is 
less than ideal, with 00 Ser in one corner of the field, and EC 53 in the bad quadrant of 
the NIRC2 detector, where one of the 8 output amplifiers of the Aladdin detector array is 
not operable and several others shows excess noise. The images were processed using various 
IRAF routines (see Section 2.2 for a general description of the data reduction procedure) 
and were astrometrically de- warped using the IDL routines provided by B. Cameron on the 
NIRC2 website 0. 

We have obtained proper motion measurements of the various H2 shock fronts from 
the H2 S(l) narrow-band images obtained in 2007 and 2009, because these images had the 
best signal-to-noise ratio for the faint shock fronts yet still showed a sufficient number of 
stars to use as an astrometric reference. By comparison, the K-hand images had more S/N 
on the stars, but the shock fronts were generally too faint to get good astrometric results. 
Establishing an astrometric reference grid on a small image near the center of a nearby 
molecular cloud is fraught with difficulties and has to rely on some simplifying assumptions. 
There are no stars in this image (Fig. 1) that could be confidently identified as distant 
background stars with negligible proper motion. The few stars in the image could possibly 
be highly reddened background stars, but they could also be faint, partly embedded young 
stars within the Serpens NW star-forming region. The best that can be done in this situation 
is to use all available stars as an astrometric reference. The images at the two epochs 
were astrometrically matched using the IRAF tasks geomap and geotran. Since the images 
had already been de-warped, we restricted the geomap solution to linear terms, to avoid 
instabilities due to the small number of reference stars. Proper motions were then measured 
using the IRAF task xregister that computes the cross-correlation between the images at the 
two epochs in a sub-image containing a feature (star or shock front) of interest. 

We summarize our proper motion measurements in Fig. 1, where the measurements 
of stars are represented by amber (light grey) arrows, while those of H2 shock fronts are 
indicated by blue (dark grey) arrows. All the stars have small residual proper motions against 
the reference system established by the ensemble of these same stars, and these proper motion 
vectors represent the combined error of the individual proper motion measurements and the 
uncertainty of the astrometric reference system defined by that ensemble of stars, as well as 
possibly genuine individual proper motion. 



http://www2.keck.hawaii.edu/inst/nirc2/forReDoc/postjobserving/dewarp/ 
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2.2. Photometric Monitoring 



The study of the temporal evolution of the reflection nebula associated with 00 Ser 
and the variability of EC 53 made use of older observa tions obtained at the UH 2.2 m 



telescope with the QUIRC camera (IHodapp et a. 



in the earlier papers on the discovery of 00 Ser (IHodapp et al.lll996al) and the initial study 



1996b[) that had al readv been published 



of the 00 Ser and EC 53 light curves jHodapplll999h . To these light curves, we have now 



added newer QUIRC data obtained up to 2004 when that car nera was decommission ed, data 
obtained at the UH 2.2 m telescope with the SIRIUS camera (INakajima et al.ll2002l ) in 2000, 
as weU as data obtained at UKIRT with WFCAM in 2008 and at CFHT with WIRCAM 
in 2011. We have also included new data fror n the recently c ommi ssioned 0.8 m Infra-Red 
Imaging Survey (IRIS) telescope described by IHodapp et al.l (120 lOl ). Finally, we have used 
the oldest available near-infrared image of the Serpens NW region that we had o btained on 
1991, Aug. 5 with the UH NIC MOS-3 camera JHodapp. Ravner. k Irwin I ll992[) at the UH 
2.2 m telescope in the K filter (jWainscoat &: Cowie Ill992l ). We used the same photometric 
apertures as for the other data sets and the same reference stars. The correction between K 
and K was determined from near-simultaneous images obtained in those two filters during 
the outburst of 00 Ser in 1995. The lowest spatial resolution data, the 1991 UH NICM0S3 
camera image and the 2010 IRIS 0.8 m telescope data, were not used for the sub-regions in 
the reflection nebula. The seeing-limited observations are listed in Table 2. 

For the study of the photometric behavior of individual knots of nebulosity within the 
00 Ser bipolar nebula, we have also included the two Keck iT-band images (Table 1), after 
scaling them to the common astrometric solution of the seeing-limited images. The Keck 
images could not be used to measure the integral brightness of the reflection nebula, because 
their field of view did not contain the integration box for the integral brightness completely. 

The images were flatfield corrected using differential dome flats, bad pixel masks were 
derived from the dome flats. The dithered object images were median combined into sky 
frames and subtracted from the flatfielded images. The resulting images were aligned and 
co-added while eliminating bad pixels in the process. After this standard data reduction, 
these seeing-limited images were transformed to a common astrometric solution using the 
IRAF tasks geomap and geotran. The available data were diverse in the total exposure 
time achieved, prevailing seeing, and photometric quality of the nights. For most, sepa- 
rate standard star observations were not available. Therefore, the photometric calibration 
relied on an ensemble of field stars. Of these field stars, we selected the four stars that 
showed the smallest scatter of the photometric measurements in the 1994 to 1998 time 
interval as secondary standar ds and calibrated them against brighter stars in the 2MASS 
catalog (ISkrutskie et al.l 120061 ). Obviously, our procedure makes the assumption that these 



selected stars were constant throughout the two decades for which photometry is being dis- 
cussed here. We have not tried to correct for color terms between the filter of 2MASS 



(ICohen. Wheaton. fc Megeath 



20031) . the classi cal ( iJohnson Ill966l ) and the Mauna-Kea- 



standard K (iTokunaga. Simons, fc Vacca II2002I ) filters used in the UH QUIRC camera and 
the other cameras occasionally used, because such effects are much smaller than the variabil- 
ity amplitudes discussed here. To avoid problems caused by the ubiquitous extended emis- 
sion, photometry on rectangular sub regions of the bipolar nebula associated with 00 Ser 
were obtained with the IRAF task imstat. The sky background was separately measured in 
similar rectangular regions selected to be unaffected by stars or extended flux. 

In a diverse data set combining data taken with different telescopes and cameras, dif- 
ferent exposure times and under varying conditions, the internal error of a photometric 
measurement underestimates the uncertainty of the result. In our case, for the higher qual- 
ity data obtained with QUIRC at the UH 2.2 m telescope, internal errors are only of the 
order of 10~^ Jy. For the UH 2.2 m telescope QUIRC and SIRIUS data, the rms scatter 
of the individual secondary standard star flux values over time is ~ 3 x 10~^ Jy. The IRIS 
0.8 m telescope data have higher rms scatter of ~ 5 x 10~^ Jy, due to poorer sampling and 
smaller signal-to-noise ratio. With these total errors, the variations in the light curve for 
different regions of the 00 Ser reflection nebula are highly significant. 



3. Results and Discussion 



3.1. Proper Motion of Ha Shock Fronts Near EC 37 = V370 Ser 



Figure 1 shows the proper motion vectors measured from the Keck NIRC2 H2 S(l) 
adaptive optics images taken in 2007 and 2009. For many of the clearly defined bow shock 
systems, these results essentially confirm the earlier results by iHodappI (119991 ) that were 
based on seeing-limited images over a 4 year interval. The new data with their superb 
spatial resolution and depth enable, for the first time, a study of the much fainter emission 
in the general a rea of EC 37, which is listed as MH02218 in the catalog of H2 shock fronts 
by lDavis et al. fl2010l)R Figure 2 is a version of the 2009 Keck AO image where an intensity- 
scaled version of the K image was subtracted, to best bring out the H2 shock fronts by 
subtracting out the continuum sources. 

The young stellar ob ject EC 37 has an average H-K f»3.6 from the photometry reported 
by lHodapp et al.l ( 1l996al ). making it one of the reddest of the point sources in the catalog of 



^The catalog is available at 



http: / /www .jach.hawaii.edu/UKIRT /MHCat 
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Serpens infrared sources by lEiroa fc Casalil ( 119921 ) . As already s howri by iHodapd (Il999l ) it 
exhibits the irregular variability characteristics of T Tauri stars fjJoy Ill945l ). In addition, it 
coincides with a s mall u nresolved knot of 3 mm emission at flux level 3.6 mJy in the map of 
Testi fc Sargent I (Il998l ). indicating the presence of substantial amounts of cold dust around 
the star, and therefore a young age of the object. On the adaptive optics H2 1-0 S(l) images, 
individual knots of emission are visible in the immediate vicinity of EC 37, primarily to the 
south of the star, in addition to more smoothly distributed flux. Very faint, knotty emission 
can be seen up to about 5" to the south of EC 37. Outside of this area, the brightest H2 
shocks are found to the east and north of EC 37. 

The large system of H2 S(l) shock fronts near EC 37 (MHO 2218) shows a pattern of 
proper motion vectors that can best be understood as an outflow emerging from EC 37. 
Th e emissio i i kno ts east of EC 37 were already imaged and their proper motion measured 
by iHodappI (119991 ). For those brighter knots, the older measurements agree very well with 
the new adaptive-optics-based proper motion vectors. In addition, fainter emission knots to 
the south-east of EC 37 could now be measured that show unambiguously that those knots 
move away from EC 37 and not from 00 Ser, which is located even further south. A few of 
the small emission knots in the immediate vicinity of the unresolved star EC 37 also have 
measured proper motion now, showing a pattern of motion radially away from the star. 

This system of emission knots near EC 37 is clearly not seen in an orientation close to the 
edge-on orientation that is favorable for seeing the full extent of a jet. The comparatively low 
extinction towards the central star that makes it directly visible at near-infrared wavelengths, 
the fact that H2 S(l) emission knots are seen all around EC 37, and the complex pattern 
of proper motion vectors all support the notion that this outflow is seen nearly along the 
direction of its axis. This is also supported by the general lack of morphological features that 
would be indicative of an edge-on orientation, such as a cometary or bipolar morphology. In 
Serpens NW, both EC 53 and 00 Ser are examples of the latter morphology. 



3.2. GO Ser 

The deeply embedded eruptive variable source 00 Ser shows nebulosity of bipolar mor- 
phology, with the western lobe being brighter than the eastern. At the spatial resolution of 
our adaptive optics co-added images of RilOO mas FWHM (see Table 1), the object at the 
center of the bipolar nebula (in the south-western corner of Fig. 1) appears slightly extended 
and elongated with the shape of a cometary nebula opening to the west. Even on these deep 
adaptive optics images, we are not seeing direct light from the central star of the 00 Ser 
reflection nebula. The whole reflection nebula brightened after the outburst of the central 
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object in 1995. However, the changes in brightness were not spatially uniform. Motivated by 
these tentative findings, we have re-analyzed older, seeing-limited images from the time of 
maximum brightness in 1995 to images taken in 2008. The images were spatially registered 
to the image taken near maximum brightness in July of 1995 using the IRAF tasks geomap 
and geotran. After adjusting the sky background to zero and after scahng of the images 
by the integrated fiux of the refiection nebula to account for the varying average brightness 
of the extended object, difference images were computed. These are shown in Fig. 3 and 
illustrate that the bipolar nebula did not change its brightness spatially uniformly. In Fig. 3, 
a faster than average decline in brightness is represented by dark tone while areas declining 
more slowly appear bright. Temporal variations in the fiux distribution of refiection nebulae 
are not a new phenomenon. For example, the case of R Mon has bee n studied as early as 



Hubble I (119161 ) . Examples observed in the near-infra red include Cep A (IHodapp &: Bressert . 
2OO9I ) and L483 (jConnelley. Hodapp. &: Fuller II2009I ) The interesting aspect of these varia- 



tions is that they point to significant variations of the extinction in the optical path of the 
scattered light on time scales of a few years, making the spatial scales of the obscuring dust 
comparable to those of our own Solar System. 

We have also extracted photometry of the 00 Ser nebula and its central extended 
source on the available seeing limited images, as described in section 2.2. Figure 4 shows 
the resulting light curves and illustrates the location of the rectangular extraction apertures 
used. In Fig. 4, the total brightness of the full 00 Ser refiection nebula, i.e., the largest 
extraction aperture, is plotted on a different flux scale than the smaller sub-apertures. It 
is clear from Fig. 4 that the overall brightness of the nebula reached a relative minimum in 
2002 or 2003, and has been slowly rising again after that time, albeit with substantial local 
variations. This result was already in dicated by tw o inde pendent K^-band photometric data 
points obtained in 2004 and 2006 bv lKospal etaP J2007h (their Fig. 4). 



Up to now, the earliest published p hotometric data point i n the 00 Ser light curve had 
been the 1994 QUIRC image shown in iHodapp et al.l (Il996al ). We have now re-analyzed 
the oldest sufficiently deep image of 00 Ser available to us. This 1991 image obtained at 
the UH 2.2m telescope with the UH NICM0S3 camera data gives a 5a detection of the 
integrated flux of the 00 Ser nebula in pre-outburst condition of K=16.8. This photometric 
value is about 2.5 magnitudes fainter than the value of i^=14.3 mag for the integrated flux 
measured in the same aperture on the 1994 QUIRC image. This demonstrates that by the 
time of the 1994 observations, the brightening process had already been under way. We have 
also double-checked this result qualitatively by the artiflcial data experiment of adding an 
appropriately intensity-scaled and smoothed version of an image of 00 Ser in its bright state 
to the 1991 image. This experiment showed that 00 Ser would have been easily visible in 
the 1991 image had it been at the same brightness as in 1994. Overall, this clearly shows 
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that a significant brightening of 00 Ser began before 1994 and maybe as early as 1991. 

Our images in the 1-0 S(l) hne of H2 show shock excited emission in the eastern lobe of 
the 00 Ser bipolar nebula. A thin arc of emission appears to outline the edge of the eastern 
outflow cavity of 00 Ser in our H— 2 1-0 S(l) narrow-band image in Fig. 1. The continuum- 
subtracted H— 2 image in Fig. 2 shows that this arc ends in a shock front dominated by line 
emission. Still further to the east, several individual shock fronts are visible whose proper 
motion (Fig. 1) traces them back to 00 Ser. This proves that 00 Ser is an accreting source 
in the early evolutionary state that is characterized by shock-excited H2 emission. The H2 
shock fronts associated with 00 Se r have now been h sted as MHO 3245 in the catalog of 



molecular hydrogen shock fronts by iDavis et al.l (I2OIOI ) 



The variable reflection nebulosity extends to the west up to the apparent edge of the 
dense molecular cloud, as judged from the density of background stars. To the east, the vari- 
able reflection nebulosity extends a little beyond the fine hne of S(l) emission just discussed, 
but is, overall, much less extended than the reflection nebulosity in the western direction. 
This morphology strongly suggests that the western lobe, i.e., the brighter reflection nebula, 
is slightly oriented towards the observer, while the eastern lobe with the H2 shock fronts 
points away from the observer and into denser parts of the Serpens NW molecular cloud so 
that reflection nebulosity becomes un-observable. 



3.3. EC 53 = V371 Ser 



The cometary nebula EC 53 shows many of the characteristics of a Class I young 
sta r seen in a nea rly edge-on orientation of its accretion disk. This object has been found 
by iHodappI (119991 ) to be variable, with strong indication of periodicity. We used the Period 
Analysis Software (Peranso) package written by T. Vanmunster, and several of the algorithms 
implemented there for a period analysis of all our seeing-limited observations of V371 Ser 
in the past 20 years. The dor ainant p e riod f ound by all algorithms is 543 days, close to 
the period initially reported in iHodappI (Il999l ). All the observations are plotted in a phase 
diagram with the 543 day period solution and shown in Fig. 5. Additional observations are 
scheduled with the IRIS system to better define the shape of the light curve. 

Our K-h&nd adaptive optics images taken with Keck/NIRC2 show this object in two 
very different phases of its light curve. In 2007, the object was near minimum while in 2009, 
it was near maximum brightness. The top panel of Fig. 6 shows the S(l) narrow-band 
image in 2007, near minimum brightness. The middle panel is a difference image of the 
S(l) line image minus a scaled K broad-band image, to bring out the H2 jet near the center 
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of the reflection nebula. The bottom panel is a scaled difference image of the 2007 - 2009 
data, showing the regions varying in phase with the illuminating source in bright tones, while 
constant (or un-phased) regions are shown in dark tone. The effectively constant parts of the 
image include the companion object (EC 53 B), and the line emission of the H2 jet. Further 
along the axis of this jet, ou tside of the field of our images present ed here, a terminal shock 
front was first identified bv iHerbst. Beckwith. fc Robbertd (119971 ) as their object Sll, and 
is included in the MHO catalog bv lOavis et al.l fcoioh as MHO 2221. In the 2007 images, 
both in broadband K and in the narrow S(l) filter, it is clearly seen that that the central 
star in EC 53 is a binary. The brighter component of this binary system (EC 53 A) is the 
variable. In its brighter state, which was observed in 2009, it so dominates the total fiux 
that it makes the detection of the companion star (EC 53 B) very difficult. A similar study 
of the uniformity of brightness changes to that done on 00 Ser was done, but essentially 
resulted in the detection of the nearly constant S(l) line emission component of the extended 
fiux. Apart from this, the change in brightness appeared uniform throughout the refiection 
nebula. 

EC 53 must be a physical binary. The companion star is located at a projected distance 
of 296 mas (92 AU) from the brighter component and has an apparent magnitude of K=19.0. 
The nearest field star to EC 53 detected in our images is 14" away. With such a low 
density of unrelated stars, the probabili ty that a 0.3" c o mpan ion is a chance alignment is ~ 
1/ 2000. From the meas urement of lEiroa &: Casalil (Il992l ). and with the extinction law 
of iRieke fc Lebofskv I (119851 ). this gives an absolute /T- magnitude of 9.75 for the companion. 
The extinction value is based on the near-infrared color excess relative to normal stars of 
the integrated light of the EC 53 nebula. The extinction towards the companion star may 
be significantly different, so the following discussion should be taken with some caution. 



The Baraffe COND03 models (IBaraffe et al.ll2003l ) for a 1 Myr star (the earliest models) 
give Mk = 9.73 for a 0.005 Mq object, while the DUSTY models give a mass between 0.004 
and 0.005. The age of EC 53 is not well determined, but based on its association with 
strong "cometary" nebulosity and a molecular hydrogen jet, and the age of its immediate 
neighbors in Serpens NW, it is probably younger than 1 Myr. The uncertainties of the 
ev olutionary models ar e substantial, in particular at the very youngest ages, as pointed out 
by iBaraffe et al.l (|2002| ). We therefore take this mass as a very rough estimate and conclude 
that the companion must be a future sub-stellar object, quite likely of planetary mass. At 
a projected distance of 92 AU, this low-r nass companion is t oo far away from the main star 
to have been formed by core accretion (IMarley et al.l 120071 ) at its present distance. Also, 
less than 1 Myr is not enough time to form giant planets by this mechanism. We therefore 
assume that the companion has formed at the same time as the main component by disk or 
core fragmentation. 
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In their study of optical photometry of T Tauri stars, iHerbst. Herbst fc Grossman 



( 119941 ) hst three main causes for variabihty in these stars: rotational modulation of dark 
star spots, modulation of bright accretion hot spots, and variable extinction due to eclipses 

by disks. While variability is common in young stars, periodicity is mostly observed on the 

timesc ale o f days to a few weeks associated with stellar rotation. For example, ICarpenter. Hillenbrand, fc 
( I2OOII ) and iHerbst et al.l ( 12002| ) have studied variability in the Orion Nebula at infrared and 
optical wavelengths, respectively. While a large fraction of stars in their respective sam- 
ples showed variability, the most common type of variability was rotationally modulated 
variations of spotted T Tauri stars with typical periods of order of several days to about 
a week, and small peak-to-peak amplitudes of a fraction of a magnitude. The short time- 
span of most variability monitoring campaigns biases the results against finding long-period 
variables. Even after considering this bias, longer period periodic variations in young stars 
appear to be very rare. One e xample of longer period variability is WL4 in p Oph, discov- 
ered by iPlavchan et al.l ( l2008l ) in the 2MASS Calibration Point Source Working Database. 
This star has a f»0.4 mag amplitude in Kg and a period of 130.87 da ys caused by eclipses o f 
each of two binary components by a circum-binary disk. For WL4, iPlavchan et al.l (120081 ) 
specif ically argue against a periodically driven accretion model for this star. On the other 
hand. iMathieu et al.l (119971 1 interpreted the photometric variations of DQ Tau (period 15.80 
days) as evidence for "varia ble accretion regulate d by the binary orbit" . Such variable accre- 
tion was also postulated by iBouvier et al.l ( 120031 ) as part of the explanation of the variations 
(period 8.2 days) observed in AA Tau. Among all the YSOs with established light curves, 
the long period, substantial amplitude and association with ongoing accretion makes EC 53 
= V371 Ser quite a unique object. 

Given that EC 53 is an outflow source, we expect that most of its luminosity is de- 
rived from accretion luminosity, and not from a, possibly pulsating, stellar atmosphere. We 
therefore believe that the most likely explanation for the periodic variations of EC 53 are 
periodic variations in the disk accretion rate, modulated by the presence of a clos e compan 



i on st ar on an eccentric orbit. This mechanism has been discussed in detail by iReipurth 



(I2OOOI ) as an explanation for FUor outbursts and for giant, well-collimated outflows. In this 
scenario, an originally non-hierarchical multiple system evolves into a hierarchical system of 
a tightly bound pair and a third component at a much larger distance that may ultimately 
be come unbound. This m odel has recently been substantiated by detailed orbital simulation 
by iReipurth et al.l (120101 ) . In a tight binary system, the main reservoir of material available 
for accretion lies outside of the binary orbit, explaining the long-wav elength spectral energy 



distrib ution peak for young accreting binaries. It has been shown by lArtymowicz fc Lubow 



(Il9961 ) that streams of material can flow through the binary orbits into the immediate vicin- 
ity of either of the two close binary components, from where they can be accreted onto the 



- 12 - 



stars by maKne tospheric accretion. These models have been detailed out by recent numerical 
simulations by Ide Val-Borro et al.l ( 1201 ll ). Even though the relationship between binary or- 
bit and accretion streams is not trivial, the model calculations show that in eccentric orbits, 
the arrival of the accretion streams at the star roughly coincides with the periastron phase 
of the binary orbit. In our case of EC 53, the observed companion (EC 53 B) at 296 mas 
projected separation is the distant third component, while the two components responsible 
for the 543 day period of light variation form the closely bound system EC 53 A. The forma- 
tion of a tightly bound binary may lead to the release of the third component of the original 
non-hierarchical system, in the case of EC 53 a free-floating object of maybe only planetary 
mass. 



So far, no repetitive FUor events have been found. Some of the s maller EXor outbur sts 
have been f ound to be repetiti ve, for example the prototypical EX Lupi (lAspin et al.ll2010l ) or 



V1647 Ori (lAspin et al. 



20061 ). but there are no cases with a clearly established periodicity. 



EC 53 may therefore represent the short-period extension of the FUor and EXor phenomenon. 



4. CONCLUSIONS 

We have presented the results of Keck/NIRC2 laser adaptive optics imaging of three 
young, accreting objects in the Serpens NW region. 

Much of the shock-excited H2 emission between the eruptive variable 00 Ser and the 
irregular variable EC 37 = V370 Ser originates in the latter. The proper motion pattern and 
the relatively low extinction to this object suggest that the line of sight towards this object 
is closely aligned with the outflow axis. 

We have identifled small patches of H2 emission east of 00 Ser, associated with the 
bipolar nebula centered on this deeply embedded outburst star. These shock fronts demon- 
strate that 00 Ser has ongoing outflow activity and is therefore still accreting mass. Based 
on images spread over many years, we show variations in the light distribution of the western 
lobe of the 00 Ser bipolar nebula that are similar to other cases of "shadow play" in young 
reflection nebulae, indicative of spatially and temporarily variable extinction near the star 
in the line of sight of the reflected light. After the 1995 outburst of 00 Ser, both the central 
source and the integrated light of the surrounding bipolar reflection nebula faded until about 
2003. After this time, the overall decline has stopped. Within the range of the spatial vari- 
ations in brightness distribution observed before, some areas of the reflection nebula have 
brightened, while others have continued to fade. 

The variable cometary nebula EC 53 = V371 Ser has an H2 jet, indicating ongoing 
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outflow activity. The earlier suggestion that the brightness variations of this object are 
periodic has been confirmed by many additional observations and a period of 543 d is now 
well established. V371 Ser has a faint, non-variable companion at a projected distance of 
92 AU that is probably a very low mass future brown dwarf or even an object of only 
planet ary mass. The v ariable cometary nebula EC 53 is being interpreted in the framework 
of the iReipurth I ( 120001 ) model whereby the ejection of a third component from a previously 
non-hierarchical system leads to the formation of a tight binary on an initially eccentric orbit 
that then imposes its orbital period on the on-going disk accretion, leading to a periodicity 
in the accretion luminosity. 



Some of the data presented herein were obtained at the W.M. Keck Observatory, which 
is operated as a scientific partnership among the California Institute of Technology, the 
University of California and NASA. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. 

The 2.2m telescope on Mauna Kea is operated by the University of Hawaii. 

The United Kingdom Infrared Telescope is operated by the Joint Astronomy Centre on 
behalf of the Science and Technology Facilities Council of the U.K. 

This paper is based on observations obtained with WIRCAM, a joint project of CFHT, 
Taiwan, Korea, Canada, France, at the Canada France-Hawaii Telescope (CFHT), which is 
operated by the National Research Council (NRC) of Canada, the Institute National des 
Sciences de I'Univers of the Centre National de la Recherche Scientifique of France, and the 
University of Hawaii. 

This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center/ California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 

The Infrared Imaging Survey (IRIS) is a joint project of the University of Hawaii, the 
Ruhr Universitat in Bochum, Germany, and the Universidad Catolica del Norte in Antofa- 
gasta, Chile. We thank L.-S. Buda, M. Dorr, H. Drass, and V. Hoffmeister for their help with 
operating the IRIS telescope. We thank B. Reipurth for helpful discussions and a careful 
reading of the manuscript. We thank the referee, Chris J. Davis, for helpful comments. 

Development of the IRIS camera was supported by the National Science Foundation 
under grant AST 0704954. 



- 14 - 



REFERENCES 

Adams, F. C, Lada, C. J., & Shu, F. H. 1987, ApJ, 312, 788 
Artymowicz, P. & Lubow, S. H. 1996, ApJ, 467, L77 

Aspin, C, Barbieri, C, Boschi, F., Di Mille, F., Rampazzi, F, Reipurth, B, & Tsvetkov, M. 
2006, A&A, 132, 1298 

Aspin, C, Reipurth, B., Herczeg, G. J., Capak, R 2010, ApJ, 719, 50 

Baraffe, I., Chabrier, G., Barman, T. S., AUard, F., & Hauschildt, R H. 2003, A&A, 402, 
701 

Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 2002, A&A, 382, 563 

Bouvier, J., Grankin, K. N., Alencar, S. H. P., Dougados, C., Fernandez, M., Basri, G., 
Batalha, C., E. Guenther, Ibrahimov, M. A., Magakian, T. Y., Melnikov, S. Y., 
Petrov, P. P., Rud, M. V., & Zapatero Osorio 2003, A&A, 409, 169 

Carpenter, J. M., Hillenbrand, L. A., & Skrutskie, M. F. 2001, AJ, 121, 3160 

Cohen, M., Whcaton, Wm. A., & Mcgcath, S. T. 2003, AJ, 126, 1090 

Connelley, M. S., Hodapp, K. W., & Fuller, G. A. 2009, A J, 237, 3494 

Davis, C. J., Cell, R., Khanzadyan, T., Smith, M. D., & Jenness, T. 2010, A&A, 511, 24 

de Lara, E., Chavarria-K., C, & Lopez-Molina, G. 1991, A&A, 243, 139 

de Val-Borro, M., Gahm, G. F., Stempels, H. C, & Peplinski, A. 2011, MNRAS, 413, 2679 

Eiroa, C, & Casah, M. M. 1992, A&A, 262, 468 

Evans, N. J., et al. 2009, ApJS, 181, 321 

Herbig, G. H. 1977, ApJ, 217, 693 

Herbst, T. M., Beckwith, S. V. W., & Robberto, M. 1997, ApJ, 486, L59 
Herbst, W., Herbst, D. K., & Grossman, E. J. 1994, AJ, 108, 1906 

Herbst, W., Bailer- Jones, C. A. L., Mundt, R., Meisenheimer, K., & Wackermann, R. 2002, 
A&A, 396, 513 

Hodapp, K. W., Rayner, J. , & Irwin, E. 1992, PASP, 104, 441 



- 15 - 



Hodapp, K.-W., Hora, J. L., Rayner, J. T., Pickles, A. J. & Ladd, E. F. 1996, ApJ, 468, 861 

Hodapp, K.-W. et al. 1996, New Astronomy 1, 177 

Hodapp, K.-W. 1999, AJ, 118, 1338 

Hodapp, K. W., & Bressert, E. 2009, AJ, 137, 3501 

Hodapp, K. W., Chilli, R., Reipurth, B., Murphy, M., Lemke, R., Watermann, R., Jacobson, 
S., Bischoff, K., Chonis, T., Dement, K., Terrien, R., & Provence, S. 2010, Proc. SPIE 
7735-45. 

Hubble, E. 1916, ApJ, 44, 190 

Johnson, H. L. 1966, ARAA, 4, 193 

Joy, A. H. 1945, ApJ, 102, 168 

Kospal, A., Abraham, P., Prusti, T., Acosta-Pulido, J., Hony, S., Moor, A., & Siebenmorgen, 
R. 2007, A&A, 470, 211 

Marley, M. S., Fortney, J. J., Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2007, ApJ, 
655, 541 

Mathieu, R. D., Stassun, K., Basri, G., Jensen, E. L. N., Johns-KruU, C. M., Valenti, J. A., 
& Hartmann, L. W. 1997, AJ, 113, 1841 

Nakajima, Y., Tamura, M., Nagashima, C, Nagayama, T., Baba, D., Sugitani, K., Nakaya, 
H., Nagata, T., Kato, D., Kurita, M., Sato, S., Jiang, Z., Ita, Y., Tanabe, T., Nakata, 
Y., Naoi, T., Oasa, Y. 2002, Proc. SPIE 4836, 29-34 

Plavchan, P., Gee, A. H., Stapelfeldt, K., & Becker, A. 2008, ApJ, 684, L37 

Reipurth, B. 2000, AJ, 120, 3177 

Reipurth, B., Mikkola, S., Gonnelley M., & Valtonen, M. 2010, ApJ, 725, L56 
Rieke, G. H. & Lebofky, M. J. 1985, ApJ, 288, 618 

Skrutskie, M. F., Cutri, R. M., Stiening, R., Weinberg, M. D., Schneider, S., Carpenter, J. 
M., Beichman, C, Capps, R., Chester, T., Elias, J., Huchra, J., Licbcrt, J., Lonsdale, 
C., Monet, D. G., Price, S., Seitzer, P., Jarrett, T., Kirkpatrick, J. D., Gizis, J., 
Howard, E., Evans, T., Fowler, J., Fullmer, L., Hurt, R., Light, R., Kopan, E. L., 
Marsh, K. A., McCallon, H. L., Tam, R., Van Dyk, S., & Wheelock, S. 2006, AJ, 131, 
1163 



- 16 - 



Testi, L., & Sargent, A. I. 1998, ApJ, 508, L91 

Tokunaga, A. T., Simons, D. A., & Vacca, W. D. 2001, PASP, 114, 180 
Wainscoat, R. J. & Cowie, L. L. 1992, AJ, 103, 332 

Wizonowich, P. L., Le Mignant, D., Bouchez, A. H., Campbell, R. D., Chin, J. C. Y., Contos, 
A. R., Van Dam, M. A., Hartman, S. K., Johansson, E. M., Lafon, R. E., Lewis, H., 
Stomski, P. J., & Summers, D. M. 2006, PASP, 118, 297 



This preprint was prepared with the A AS macros v5.2. 



-17- 



Table 1. Log of Keck AO Observations 



Date (UT) 


JD 


Filter 


t 


nreads 


coadds 


frames used 


FWHM [mas] 


2007-06-10 


2454262 


H2 


30s 


8 


1 


125 


88 


2007-07-09 


2454291 


K 


2s 


8 


45 


25 


123 


2009-06-20 


2455003 


H2 


60s 


16 


1 


82 


114 


2009-06-29 


2455012 


K 


2s 


8 


45 


15 


97 



Table 2. Log of Observations 



Date JD Instrument Filter 



5 Aug 1991 


2448474 


UH NICMOS3 


K' 


16 Aug 1994 


2449581 


UH QUIRC 


K 


13 Jul 1995 


2449912 


UH QUIRC 


K 


12 Oct 1995 


2450003 


UH QUIRC 


K 


6 Feb 1996 


2450120 


UH QUIRC 


K 


28 Jul 1996 


2450293 


UH QUIRC 


K 


3 Oct 1996 


2450360 


UH QUIRC 


K 


25 Apr 1997 


2450564 


UH QUIRC 


K 


17 Aug 1997 


2450678 


UH QUIRC 


K 


5 Feb 1998 


2450850 


UH QUIRC 


K 


30 Sep 1998 


2451087 


UH QUIRC 


K 


18 Aug 2000 


2451775 


UH SIRIUS 


K, 


1 Aug 2002 


2452488 


UH QUIRC 


Kmko 


10 Jun 2003 


2452801 


UH QUIRC 


Kmko 


31 Jul 2004 


2453218 


UH QUIRC 


Kmko 


21 May 2008 


2454608 


UKIRT WFCAM 


Kmko 


17 Aug 2010 


2455426 


IRIS 


K, 


28 Aug 2010 


2455436 


IRIS 




14 Sep 2010 


2455454 


IRIS 




1 Oct 2010 


2455471 


IRIS 


K« 


20 Oct 2010 


2455490 


IRIS 


K, 


20 Mar 2011 


2455641 


IRIS 
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Fig. 1. — Proper motion vectors measured on the Keck/NIRC2 2.12 //m images over the 
time interval 2007 - 2009. The underlying image is the 2009 Keck NIRC2 LGSAO 2.12 /im 
image. Amber (light grey) arrows indicate the proper motion of stars in the field, and are 
largely an indication of the uncertainties of the astrometric reference that is based on the 
ensemble of these stars. Blue (dark grey) arrows indicate the measured proper motions of 
H2 shock fronts. 
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Fig. 2. — Continuum subtracted v=l-0 S(l) image of the shock fronts associated with EC 37 
and 00 Ser based on the 2009 Keck data. The broad-band (continuum) image was scaled to 
subtract out the flux of the 00 Ser reflection nebula. The reddest object in the frame, the 
central condensed area of this reflection nebula, therefore appears white (over-subtracted) 
while all other, less reddened continuum objects, appear black. 
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Fig. 3. — Scaled difference images based on seeing-limited if -band images of the 00 Ser 
reflection nebulosity between 1996 and 2008. The field of view is approximately 23" x 19". 
All images used were first scaled by the integrated signal of the reflection nebula to account 
for the overall brightness variations of the object. For each of the months and years indicated 
(for prccicc dates see Table 2), the difference image is the fiux-scaled 1995 peak brightness 
image minus the image taken at the given date. In this scaling, white tone indicates slower 
fading than average, while dark tone indicates faster than average fading. 
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Fig. 4. — Light curve of the integrated K-hand hght from 00 Ser and from individual 
regions of the reflection nebulosity. This photometry is based on all the observations listed 
in Tables 1 and 2. The top (green) light curve represents the integral flux over the whole 
reflection nebula. The flux calibration for this curve is to the right side of the figure. The 
other light curves are for the smaller sub-regions of the reflection nebula, and they refer to 
the flux scale on the left, which is a factor of 5 smaller than the right scale. The lower panel 
of the figure shows, on the same flux scale as for the sub-apertures, the measurements of the 
four nominally constant reference stars, as an indication of typical total measurement errors. 
The horizontal axis gives the time both in abbreviated Julian dates and calendar years. 
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Fig. 5. — Phase diagram of the if -band variations of EC 53. The hght curve shows a rapid 
rise and slower dechne from the maximum. 
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Fig. 6.— Top panel: Keck NIRC2 LGSAO image in the 2.12/im S(l) filter of V371 Ser = 
EC 53 taken near minimal light in 2007. Middle panel: Difference image of the S(l) line 
image and a K broad-band image. The broad band image was scaled to that the reddest 
object in the field, the variable source V371 Ser A, appears white, the companion V371 Ser B 
is subtracted away, and regions with strong line emission appear dark. The line-dominated 
flux extends to the south-east from V371 Ser A, while the overall reflection nebula (top 
panel) extends straight to the south. Lower panel: Difference of 2.12//m S(l) line images 
taken in 2007 (near minimum) and 2009 (close to maximum) of V371 Ser. The image was 
scaled to show regions varying strongly and in phase with V371 Ser in white, while regions 
of constant ffux are dark. The companion V371 Ser B is visible as a black dot, indicating 



